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ABSTRACT 

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with 
resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to 
~ 15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September 
to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science 
Verification (SV) targets. This paper presents an overview of the campaign and its main results, including 
an investigation of the short-term coherence properties and systematic phase errors over the long baselines at 
the ALMA site, a summary of the SV targets and observations, and recommendations for science observing 
strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared 
to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive 
program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as 
fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens 
up new parameter space for submm astronomy. 

Subject headings: instrumentation: interferometers—submillimeter: general—telescopes—techniques: high 
angular resolution—techniques: interferometric 
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1. INTRODUCTION 

The Atacama Large Millimeter/submillimeter Array 
(ALMA) is a millimeter/submillimeter (mm/submm) interfer¬ 
ometer located in the Atacama desert of northern Chile at an 
elevation of about 5000 m above sea level. The high-altitude, 
dry site provides excellent atmospher ic transmission over th e 
frequency range 85 GHz to 900 GHz dMatsushita et alJl999l) . 
ALMA is currently in its third year of science operations and 
was formally inaugurated in 2013 March. Until now, science 
observations have used configurations with baselines from 
100 m to ~1.5 km, with some limited testing of a ~3-k m 
baseline in 2013 (lAsaki et al.1(20141 : iMatsushita et al.ll2014l) . 
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Fig. 1.— Example LBC array configuration (in this case the array that was 
used for the 3C138 Band 6 observations in Appendix A). The black points 
show the nominal LBC antennas. The five antennas near the center (red 
points) are not part of the nominal LBC array, but were useful for measuring 
more extended emission (the number of these antennas varied; see Section 2 
for details). The axis units are in meters. 


To test the highest angular resolution capability of ALMA 
using baseline lengths of up to ~15 km at selected frequen¬ 
cies, the three-month period from 2014 September to Novem¬ 
ber was dedicated to carrying out the 2014 ALMA Long Base¬ 
line Campaign (LBC Fl The approximate resolutions that can 
be achieved with the longest baselines are 60 mas at 100 GHz, 
25 mas at 250 GHz and 17 mas at 350 GHz (but these can vary 
by ~20% depending on the imaging parameters). The major 
goal of the campaign was to develop the technical capabilities 
and procedures needed in order to offer ALMA long baseline 
array configurations for future science observations. 

This paper presents an overview of the ALMA LBC, focus¬ 
ing on the technical issues affecting submm interferometry on 
baselines longer than a few kilometers. In ij2] we describe 
the LBC array and campaign test strategy. §[3] describes the 
effects of short-term phase variation due to the atmosphere 
and a method for determining if conditions are sufficiently 
stable for imaging. In $4] we discuss the systematic phase er¬ 
rors found between the calibrator and science target. In ^5j 
an overview of Science Verification (SV) at long baselines is 
given. Images and initial science results on the SV targets 
are presented in three accompanying papers (ALMA Partner¬ 
ship et al. 2015a,b,c). An illustration of the quality of the 
ALMA calibration and imaging is given by a comparison of 
preliminary ALMA SV and Very Large Array (VLA) images 
of 3C138 with the same resolution (Appendix A). In )j6j we 
present conclusions drawn from the LBC and recommenda¬ 
tions for science observing using long baselines with ALMA. 


2. LONG BASELINE CAMPAIGN OVERVIEW 
2.1. The LBC Array 

Since many of the distant antenna pads had not been pre¬ 
viously powered or occupied, a coordinated effort was made 
from April to August 2014 to prepare a sufficient number of 
antenna stations beyond 2 km from the array center. The con- 

43 The LBC was led by the Extension and Optimization of Capabilities 
(EOC) team, which includes members from the Joint ALMA Observatory 
(JAO) Department of Science Operations. It was a collaborative effort by an 
international team including members from the JAO, the ALMA Regional 
Centers, and the JAO expert visitor program. 
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Fig. 2.— The uv distribution for the the ~l-hr 3C138 Band 6 observations. 
The M-axis is the east-west spacing and the v-axis is the north-south spacing. 
The axis units are in kilo-lambda (kA, where 1000 kA = 1300 meters). 

figuration process began with an initial test in late August 
2014 when a single antenna was moved out to a 7 km baseline. 
The nominal LBC configuration consisted of 21-23 antennas 
on baselines of between 400 m and 15 km and was avail¬ 
able from the end of September until mid-November 2014 
(with the two longest baseline antennas being added in mid- 
October). In addition, typically 6-12 antennas were available 
on baselines less than 300 m that were useful for imaging the 
more extended sources (though since they were not part of 
the nominal LBC configuration, the number of these antennas 
on short spacings varied from day to day and with observing 
Band). Thus, the total number of antennas used during the 
campaign typically ranged from 22-36, depending on observ¬ 
ing date and observing Band. An example configuration used 
during the campaign (in this case for the SV observations of 
3C138; see Appendix A) is shown in Fig. |T] The resultant u-v 
coverage for a ~l-hr observation of 3C138 with this array is 
shown in Fig. [2] 

2.2. LBC Test Strategy 

The norm al calibration mode for AL MA observing is phase 
referencing dBeaslev & Conwavil 19951) . Over the length of an 
experiment that can last for several hours, this observing mode 
alternates short scans of the science target and a nearby quasar 
that is used to calibrate the target data. Hence, the outcome 
of the long baseline observations depends strongly on the ac¬ 
curacy with which the phase measured on the calibrator can 
be transferred to the target. The LBC concentrated on the ac¬ 
curacy of this transfer by: (1) performing test observations of 
quasars to establish the properties of the phase coherence of 
the array over long baselines; (2) determining how to opti¬ 
mize observing strategy to achieve good imaging results; and 
(3) observing, calibrating, and imaging SV targets and other 
test targets to demonstrate the end-to-end capability of ALMA 
long baseline observations. 

Key plans for the LBC testing included: (1) Source stares: 
30-min observations of a single bright source to determine 
the temporal phase variation statistics as a function of base¬ 
line length; (2) Short phase reference tests: alternating ob¬ 
servations of two close sources to determine the accuracy of 


the phase transfer and subsequent image errors; (3) Go/noGo 
tests: development of an online method to determine the 
near real -time feasibility of long baseline observations (Sec¬ 
tion EH. (4) Cycle time tests: phase referencing tests with 
different intervals between calibrator scans; (5) Baseline de¬ 
termination: observations of many quasars distributed over 
the sky for 30 to 60 min to determine antenna positions and 
delay model errors; (6) Weak calibrator survey: measur¬ 
ing the flux density of candidate calibrators for suitability as 
phase reference sources; (7) Calibrator Structures: imaging 
of calibrators at long baselines to search for significant angu¬ 
lar sizes; and (8) Astrometry: phase referencing among many 
close quasars to measure the long baseline source position ac¬ 
curacy. 

Most test observations were made at 100 GHz (ALMA 
Band 3). The observed phase fluctuations are associated with 
variations in propagation time (delays) in the ALMA sys¬ 
tem or in the atmosphere, which are also described as path- 
length variations. The propagation changes are generally 
non-dispersive so that the phase fluctuations will scale with 
frcqucnc\Fl (although there are significant dispersive effects 
in the contributions due to water vapor at some frequencies 
above 350 GHz; these effects can be estimated). 

3. SHORT-TERM COHERENCE 

Imaging using phase referencing techniques requires a rea¬ 
sonably phase-stable array. Hence, an early goal of the LBC 
was to determine the short-term (5 to 60 sec) phase rms prop¬ 
erties of ALMA over a variety of conditions. In addition to 
phase noise, systematic phase offsets between the science tar¬ 
get and calibrator were found; in (JfQ we describe their origin 
and how they were minimized. 

One of the main contributions to phase instability at mm 
wavelengths is the fluctuation of the amount of water vapor in 
the atmosphere. The ALMA site was chosen for its low aver¬ 
age water vapor content and excellent phase stability. Never¬ 
theless, at baselines longer than 1 km, the short-term phase 
variations may make imaging impossible. A good rule of 
thumb is that if the rms phase variations are a (rad), then the 
approximate loss of coherence (the decrease of the peak in¬ 
tensity of a point source caused by these random phase fluc¬ 
tuations) is exp[(-cr 2 /2)] dRichardsll2003l) . For a = 30° or 
60° the coherence is respectively 87 or 58%. Hence, a gen¬ 
eral guideline is that the loss of coherence is acceptable and 
reasonably accurate image quality can be obtained if the rms 
phase fluctuations are < 30°. 

3.1. WVR correction and the Spatial Structure Function 

To estimate the path variations associated with the water 
vapour component, each antenna is equipped with a Water Va¬ 
por Radiometer (WVR). The WVR is a multi-channel receiver 
system (Emrich et al. 2009) that makes continuous observa¬ 
tions of the emission in the wings of the 183 GHz water line 
along the line of sight to the astronomical source. A descrip¬ 
tion of this system, and of the way in which the measurements 
are used to estimate the variations in the amount of Precip- 
itable Water Vapor (PWV0 in the path to each antenna, is 

44 A useful conversion is that a path length change of 1 mm will produce 
a path delay change (assuming propagation at c) of 3.3 psec. The 1 mm path 
length change will produce a phase change of 120° at 100 GHz (Band 3), 
300° at 230 GHz (Band 6), and 420° at 340 GHz (Band 7). 

45 Each mm of PWV along the line of sight will result in a path length 
increase of 6.5mm: lThomDson et all (20011 ). 
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Fig. 3.— The spatial structure function (SSF). The phase rms 0 (square- 
root of the SSF; converted to a path length in microns) versus baseline length 
is shown for a target at three stages of reduction. The experiment was 15 
minutes in duration. The black points show the SSF for the original visibility 
data. The red points show the SSF points after applying the WVR correction 
for this source. The orange points show the SSF for this source after phase 
referencing with a calibrator that is 1.3° away from the target with a cycle 
time of 20 sec. The PWV during this experiment was 1.44 mm with a wind 
speed of 7 m/sec. 
a See footnote 44. 



Fig. 4.— Average phase rms for 5 to 15 km baselines versus coher ence. 
The rms phase was determined from a set of Go/noGo (see Section E2) 
observations that were followed by short phase referencing experiments. The 
coherence is the ratio of the phase referenced image peak density divided 
by the self-calibrated image peak flux density. The thin red line shows the 
theoretical relationship between the phase rms (radians), a, and coherence, 
i.e. exp(-cr 2 /2) for a random phase distribution (see Section[3). 

given in lNikolic et al.1 (120131) . This WVR correction typically 
removes about half the short-term phase fluctuations, and in¬ 
creases the proportion of time that phase referencing observa¬ 
tions will produce good quality images. Even in good condi¬ 
tions, however, applying a correction to the phases based on 
these estimates still leaves residual fluctuations that are much 
larger than the estimated errors (which, with clear skies and 


PWV < 2mm are believed to be less than 20 microns, although 
they can be much larger when clouds of ice or liquid water 
are present); see Figure 3. These residuals are thought to be 
mainly due to dry atmosphere (i.e. density) fluctuations (also 
see Section l4~2l i. 

The properties of the phase rms as a function of baseline 
length are important for deciding when and how to observe at 
long baselines. Fig.QJshows a typical relationship of the phase 
rms, <7, as a function of baseline length, /;, for a target at three 
stages of analysis. The cr-b relationship is called the Spa¬ 
tial Structure Function (SSF). The characteristic shape is sim¬ 
ilar for both the uncorrected data and for the WVR corrected 
data, except for the decrease of the variations by about 50%. 
For short baselines, the rms phase increase s as a ss Z> 0 83 , in¬ 
dicative of a 3-D Kolmogorov spectrum (ICarilh & Holdawavi 
Il999f) . The slope then decreases to a 2-D Kolmogorov spec¬ 
trum with dependence a ~ Z? 0,33 at about 3 km, which is 
roughly the scale height of phase turbulence. This scale height 
is an average of the wet atmosphere and dry atmosphere scale 
sizes of 1 and 5 km at the ALMA site. 

After phase referencing, the shape of the SSF is altered, as 
shown by the orange points in Fig. [3] In this example, the cal¬ 
ibrator is only 1.3° away from the target, the cycle time is 20 
sec and the integration time on the calibrator is only 6 sec. 
Only a small fraction of calibrators are sufficiently strong, 
even at Band 3, to provide adequate signal-to-noise for ac¬ 
curate phase referencing calibration in this short integration 
time. Even in the ideal case of a sufficiently strong calibra¬ 
tor, for baselines less than 1 km, there is little decrease in the 
target rms after phase referencing. However, beyond a base¬ 
line of about 1 km, the target rms becomes less dependent on 
baseline length since the phase fluctuations with scale sizes 
greater than 1 km are well correlated between the target and 
calibrator with a 20 sec switching cycle time. 

3.2. Go/noGo System 

At the beginning of the campaign, it was hoped that the 
properties of the rms phase fluctuations (both before and af¬ 
ter WVR correction) could be predicted from measurable 
weather parameters such as the average PWV, PWV rms, 
wind speed, and pressure rms. If so, then algorithms associ¬ 
ated with these measured conditions could be used to indicate 
in advance if the phase parameters are adequate for imaging at 
a specified frequency; namely, that the short-term phase rms 
would be less than about 30° for the longer baselines. This 
presumption, however, turned out to be not always true. 

A direct method to determine the current ALMA phase rms 
is from a short observation of a strong source. A simple ob¬ 
serving procedure called Go/noGo was developed, consisting 
of a 2-min observation of a strong quasar at Band 3, followed 
by online data analysis that rapidly determines the SSF with 
the WVR correction applied. To confirm that the Go/noGo 
structure function phase rms (averaged over many baselines 
between 5 to 15 km) is well correlated with phase referencing 
image quality, many Go/noGo observations that were carried 
out during the LBC were followed by short reference obser¬ 
vations of calibrator-target pairs, with a typical 3.5° separa¬ 
tion and cycle time of 60 sec. The plot of the Go/noGo rms 
phase versus image coherence from the phase referencing ex¬ 
periment is shown in Fig. [4] This demonstrates that the target 
image coherence is reasonably well correlated with the rms 
phase at the longer baselines of the calibrator. The reason for 
the somewhat lower image coherence than expected from the 
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Fig. 5.— The effect of the height difference delay term. The residual delay/phase for one baseline (after fitting for the best antenna positions) is plotted versus 
the sin(elevation) for 50 quasar scans that form a typical baseline observation. The baseline lengt h is 3 .5 km with an antenna height difference of 100 m. The red 
points show the residual delay/phases that used the nominal ALMA CALC delay model (Section l4.lj that assigned the measured pressure from the one weather 
station near the array center to both antennas. The blue points show the residual delay/phase for another baseline observation in which the estimated pressure for 
each antenna was estimated using the pressure lapse rate from the barometer near the array center. 


rms phase variations are discussed in ^4] 

4. SYSTEMATIC PHASE ERRORS 

In addition to the stochastic-like phase variations between 
the calibrator and target described in Section [3] there were 
systematic antenna-based phase offsets between the calibrator 
and science target that persisted on timescales of many min¬ 
utes to hours. These were found to be caused mostly by errors 
in the correlator delay model. The offsets were found to scale 
roughly as the calibrator-target separation, but were nearly un¬ 
affected by the cycle time. Such systematic offsets can have 
serious impact on the target image quality because they are 
persistant and produce image artifacts (e.g. large side-lobes 
and spurious faint components), in addition to the blurring of 
the target image that is associated with short-term phase fluc¬ 
tuations. 

4.1. The Delay Model 

The signals from all antennas must be combined precisely 
in phase at correlation to obtain accurate visibility phases. A 
critical part of the ALMA online control software, called the 
delay server, calculates the expected relative delay of the sig¬ 
nals between each antenna from the ALMA array parameters 
dMarson et al 12008]) . If the delay model (DM; which is calcu¬ 
lated using the CALCP1 third-party software) is accurate, the 
visibility phase for any point-like quasar with known position 
should be constant with time and independent of the quasar’s 
position in the sky. 

An important part of the DM is the estimate of the differen¬ 
tial tropospheric delay between each antenna from the source. 
As described above, the wet delay component is calculated 
from the 183 GHz emission assuming a model temperature 
profile, and is included in the DM using the WVR measure¬ 
ment. The zenith dry air delay t, above antenna i is accurately 
given by r, « 0.228 P, where P, is the dry pressure in mbars 

46 http://lacerta.gsfc.nasa.gov/mk5/help/calc_01.txt 


at the antenna ([Thompson et al.ll20011) . For an observation of 
a target at elevation e, the CALC model delay is (rf) / sin(e). 

Given that only one weather station near the array center had 
so far been available at the time of the LBCP 7 ! at ALMA, the 
estimate of the dry air delay at each antenna is not as accurate 
as desired. This inaccuracy results in antenna-based phase 
offsets that differ between a calibrator and science target and 
hence produce relatively constant phase offsets between them. 

4.2. Measurement of Delay Model Errors 

The presence of DM errors was suspected from the base¬ 
line observations that consisted of about 50 to 100 ten-second 
quasar observations distributed over the sk\PT. Many such ob¬ 
servations have been made in order to determine the accurate 
relative positions of the antennas which are frequently moved 
from one antenna pad to another as the ALMA configuration 
changes. The a priori antenna positions are usually more than 
1 mm in error, so the baseline observations provide the data 
needed to update antenna positions, generally to an accuracy 
of about 50 microns. Over a few years, it was found that the 
measured position changes of fixed antennas between base¬ 
line calibration observations, separated by several hours to a 
few weeks, were often larger than 100 microns and sometimes 
well over 1 mm for unmoved antennas that were more than 
1 km from the array center. 

These apparent antenna position changes were traced to 
the implementation of the dry ah' delay term in the CALC 
DM. Fig. 0 illustrates the results of an experiment on 2014 
September 16 with two 30 min baseline observations which 
confirmed the DM error for a 3.5 km baseline with a height 
difference of 100 m between the two antennas. One experi¬ 
ment used the DM in which the pressure at each antenna was 
set equal to that measured by the one sensor. After fitting for 

47 Installation and testing of several additional weather stations distributed 
over the array is planned for the end of 2015. 

48 http://legacy.nrao.edu/alma/memos/html-memos/alma503/memo5 



























6 


ALMA Partnership et al. 


TABLE 1 

Long Baseline Astrometric Results 


Source 

J0538-4405 a 

J0519-4546 

J0455-4615 

J0522-3627 



Sep (Deg) 

0.0 

3.8 

7.9 

8.2 



RMS (mas) b 

0.13 

0.20 

0.35 

0.08 



DATE 


RA, DEC offset (mas) 


PWV(mm) 

ELEV(deg) 

Sep 22 

-0.45, +0.15 

-0.44, -2.20 

+ 11.0,-2.96 

+0.48,+1.82 

0.7 

66 

Oct 03 

-0.09, -0.10 

+0.28, -4.57 

+4.03,-8.14 

+1.32,-2.99 

2.7 

71 

Oct 14 

+0.06, +0.08 

-2.69, +0.25 

— 

+0.74, -3.49 

0.6 

55 

Oct 14 

-0.00, -0.03 

+0.54,+1.27 

+8.81,+3.84 

-4.00, -1.70 

0.9 

71 

Nov 04 

-0.03. -0.02 

-0.21,+2.12 

+7.08, +4.80 

-4.60, +0.60 

0.9 

48 

Nov 17 

-0.01, +0.05 

-2.72,-1.24 

-6.56,-1.96 

-9.02,-0.17 

1.3 

66 

Mean 

-0.09, +0.02 

-0.87, -0.73 

-4.87, -0.88 

-2.51,-0.99 



STD (Mean) 

0.08, 0.03 

0.54, 0.92 

2.75,2.15 

1.15,0.78 




Note. — Details are given in Section l4~4l 
a Phase reference source. 

b Theoretical RMS, defined as the angular resolution divided by the theoretical SNR, where the latter is 
derived from the peak flux density of the source divided by the expected image rms noise level. 


the best baseline, the residual fit, shown by the red points, 
contains a large residual phase versus elevation term. In the 
subsequent experiments, the pressure at each antenna was es¬ 
timated using the approximate pressure lapse rate. After the 
best baseline fit to the data, the residual phase versus elevation 
is flat. Since some antennas in the long baseline array have a 
height difference from the array center of over 200 m, even 
larger systematic phase errors could be encountered. 

Without a reasonable pressure estimate for each antenna, 
the target and phase will have a systematic offset that will 
only slowly change. For example, the residual phase between 
a calibrator at elevation 55° and a target at elevation 60° is 
about 110°; this phase offset is not removed by the phase ref¬ 
erencing. After the September demonstration of the issue, the 
ALMA DM was updated to include an estimate of the pres¬ 
sure at all antennas using the lapse pressure rate and the h eight 
of the current single pressure monitor (as noted in Section FhTl 
additional pressure monitors distributed across the array will 
be available in fut ure). This height-delay co mpensation is also 
used at the VLA (iFomalont & Perlevl 19991) . 

Even after the correction of the antenna height delay dif¬ 
ferences, additional baseline observations during the last part 
of the LBC still showed apparent antenna position offsets of 
about 1-5 mm for most antennas 5 to 10 km from the cen¬ 
ter, which scaled roughly with distance from the array center. 
These apparent antenna position changes are consistent with 
the un-modeled pressure changes expected over the 15 km 
region of the Chajnantor plateau. However, by using a cali¬ 
brator close to the target this effect is minimized ; thi s requires 
a larger catalog of potential calibrators (Section |4~3T ). 

Additional observational techniques can be employed to 
model the dry term delay residuals. For example. Very 
Long Baseline Array (VLBA) observations often include a 
short baseline-type observation (20 sources in 20 min) to 
determine the residual zenith path delay over each antenna 
dReid & Honmal[i014l) . Such options may be explored for fu¬ 
ture work. 


brator catalog] in September 2014 contained 700 entries of 
quasars with positional accuracy < 2 mas from Very Long 
Baseline Interferometry (VLBI) observations and with a 100 
GHz flux density > 25 mJy. Over the ALMA sky between 
-90° and +45° declinations, the mean angular distance of an 
ALMA catalog entry from a random target is 3.5° with a 25% 
chance that the closest calibrator is > 5° away. The number 
of suitable calibrators in the catalog, especially for the long 
baseline observations, therefore needed to be substantially in¬ 
creased. To this end, a survey of weak calibrators was ini¬ 
tiated in mid-September to observe candidate sources from 
the AT20G iMassardi et all (1201 ll) and VLA calibrator cata¬ 
logs 0 to determine their flux densities at 100 GHz. This 
list of 4200 candidate sources was compiled from sources po¬ 
tentially stronger than 25 mJy at 100 GHz, and observations 
prioritized the ~3000 sources with VLBI positions^] having a 
positional accuracy of < 2 mas. Sources as faint at 10 mJy at 
Band 3 may potentially be used as phase calibrators, but find¬ 
ing the faintest acceptable calibrators will probably require 
future targeted searches around a source. 0 
About 20 of the brightest ALMA calibrators were also im¬ 
aged with the LBC array to determine if they were resolved 
at the longer baselines. Since most of the sources have been 
previously imaged using VLBI baselines of 5000 km at cm- 
wavelengths and found to be less than about 5 mas in angular 
size, it was expected that these calibrators would be nearly un¬ 
resolved sources at ALMA long baseline resolutions. Two of 
the 20 sources, however, had faint inner jets whose brightness 
was a few percent of the bright core point component, but this 
structure has little effect on their use as calibrators of ampli¬ 
tude and phase on long baselines. A few of the brighter cal¬ 
ibrators were already known to have large arcsec-scale struc¬ 
ture (J0522-3627 and 3C273); this also has no significant ef¬ 
fect on their use as long baseline calibrators. 

4.4. Astrometric Accuracy 


4.3. The Weak Calibrator Survey and Calibrator Structure 

To facilitate an optimal calibrator choice for a science tar¬ 
get, most observatories support a source catalog that contains 
information about candidate calibrators. The ALMA cali- 


49 https://science.nrao.edu/facilities/alma/aboutALMA/Technol 
http://www.eso.org/sci/publications/messenger/archive/no.155-r 

page 19 

50 http://www.aoc.nrao.edu/$\sim$gtaylor/csource.html 

51 http://astrogeo.org/vlbi/solutions/rfc/atmos_2014d 

52 https://science.nrao.edu/facilities/alma/aboutALMA/Technol 
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TABLE 2 

ALMA Long Baseline Science Verification Targets 


Target 

Coordinates 11 

Band b 

Scope 0 

N 0 „, d 

N,/ 

Lw f 

Freq. g 

Obs. Date* 1 

Id.' 

Juno 

ephemeris target 

6 

cont., ephemeris 

30-33 

5 

0.3* 

224.0, 226.0, 240.0, 242.0 

10/19 

13 

Mira 

02 ,, 19"'20L79 -02°58'39".5 

3 

SiO, cont. 

31-33 

3 

1.5 

88.2, 98.2. 100.2, 86.8, 

86.2, 85.6, 85.7 

10/17-10/25 

14 


04 ,1 31 m 38 s .45+18 O 13 , 59".0 

6 

SiO, cont. 

35-36 

3 

1.0 

229.6.214.4,214.1.215.6. 
217.1,232.7, 231.9 

10/29-11/01 

14 

HL Tau 

3 

CO, CN, cont. 

32-35 

7 

3.2 

102.9,101.1,115.3,113.5 

10/28-11/14 

15 



3 

HCN, HCO+, cont. 

33-35 

7 

3.5 

90.8,100.8,102.8,88.6,89.2 

10/14-11/13 

15 



6 

cont. 

28-36 

9 

4.7 

224.0,226.0,240.0,242.0 

10/24-10/31 

15 


05 , ’21 m 09 I .9+16°38 , 22" 

7 

cont. 

27-36 

10 

5.1 

336.5,338.4,348.5,350.5 

10/30-11/06 

15 

3C 138 

3 

cont., polarization 

27-30 

6 

2.0 

90.5, 92.5, 102.5, 104.0 

11/10-11/19 



09 ,! 03"'1P.61 +00°39 , 06".7 

6 

cont. 

29-31 

5 

1.6 

224.0, 226.0, 240.0, 242.0 

11/09-11/14 


SDP.81 

4 

CO, cont. 

22-27 

12 

5.9 

144.6,154.7,156.4,142.7 

10/21-11/11 

16 



6 

CO, H 2 O, cont. 

30-36 

9 

4.4 

228.0, 230.0, 243.0, 244.5 

10/12-11/09 

16 



7 

CO, cont. 

31-36 

11 

5.6 

282.9. 294.9, 296.9, 284.9 

10/30-11/04 

16 


NOTE. — Further details of the Juno, HL Tau, and SDP.81 observations and results are given in three accompanying papers (ALMA Partnership et al. 2015a,b,c). The data is 
publicly available from the ALMA Science PortafJ 
a Coordinates of the phase center (J2000) 

b ALMA Bands. Bands 3, 4, 6, & 7 correspond to frequencies of approximately 100 GHz, 140 GHz, 230 GHz & 340 GHz, respectively. 
c Scope and aim of the observations. These include spectral line and/or continuum imaging at high angular resolution, plus polarization and ephemeris targets. 

^ Number of antennas in the array for each execution. Typically, between 1-5 of the total number of antennas were flagged for a given execution. Note that the number and 
configuration of the antennas on very short spacings varied from day to day (see Section [2j. The number of antennas also varied with observing Band, with the fewest antennas 
available in Band 4 (due to fewer antennas with Band 4 receivers available during the LBC). 
e Total number of executions of the scheduling block. 

^ Total effective integration time on source (i.e., after flagging), in hours. *For specific details of Juno, see ALMA Partnership et al. (2015a). 

£ Mean center frequency of each spectral window (spw) in GHz. Channel widths were 15.6 MHz for continuum windows, 2.0 GHz bandwidth. Channel widths varied for spectral line 
windows. For Mira, they were 61-122 kHz, 0.059-0.117 GHz bandwidth. For HL Tau, they were 61 kHz, 0.117 or 0.243 GHz bandwidth. For SDP.81, they were 0.488-1.953 GHz, 
1.875 GHz bandwidth. 

^ Range of dates of the observations. 

1 The project code identifier of the dataset can be obtained by replacing “XX” in ADS/JAO.ALMA#2011.0.000XX.SV with the number in this column. 
a http ://w w w. almascience. org 


During the campaign, many hour-long experiments, cycling 
among three or four quasars within a radius of 10°, were car¬ 
ried out. All of the quasars have an a priori position accu¬ 
racy of < 0.3 mas, and were observed sequentially with 1-min 
scans at Band 3. Using one of the quasars as the phase ref¬ 
erence calibrator, images of the other quasars were obtained 
and the positional offset for each source was determined by 
the displacement of the quasar peak from the center of the im¬ 
age. Uniform weighting with only spacings longer than 1 km 
were used to obtain the highest resolution and most accurate 
positions. 

The results for the same quartet of quasars observed six 
times over the LBC are given in Table Q] The source J0538- 
4405 is the phase reference source, so its position should be 
close to zero. The separation of the sources from J0538-4405 
in degrees and the theoretical positional rms error in mas are 
listed in the first two rows. The subsequent rows then give 
the R.A. and deck offset for each source for each of the six 
observations, with the mean positional offset and the standard 
deviation of the mean at the bottom. The results show that 
the positional offsets of the three target sources are signifi¬ 
cantly larger than those expected from the image noise alone 
(typically ~l-5 mas). The source, J0519-4546, closest to the 
phase reference source, shows the smallest systematic offset 
(~0.8 mas). The other two sources, one to the east and one to 
the south of J0538-4405, have larger offsets. This relationship 
is consistent with that produced by the relatively systematic 
atmospheric delay model errors discussed in Section 14.21 In 
future, it will be possible to use the apparent positional off¬ 
sets of three calibrators to determine more information about 
the delay model error over the array, and then remove the er¬ 


rors to obtain more accurate positions of the calibrators. Such 
multi-calibrator observations and analyses have proved suc¬ 
cessful with the VLBA for significantly improving the as¬ 
trometric precision (iFomalont & Koneikinil20021) and are now 
being tested for ALMA. 

The nominal astrometric accuracy from the LBC tests, 
given by the average rms in Table 1 for the three sources, is 
an rms positional error of ~1.5 mas. This is for an average 
calibrator-target separation of ~ 6°, with an observing period 
of one hour, with a maximum baseline of 12 km. Given a 
sufficiently strong point source, this accuracy is independent 
of observing frequency. The predicted ALMA astrometric ac¬ 
curacy is ~0.18 mas (Lestrade 2008), assuming the use of 
WVR corrections and a typical calibrator-target separation of 
6° (which is within the range used in the LBC). However, 
this predicted value assumes that the pressure measurement at 
each antenna would be accurate to ±2 mbars. As discussed 
in Section 4, with the availability of only one weather sta¬ 
tion during the LBC, the inferred pressure for antennas many 
kilometers from the pressure sensors, using a simple plane- 
parallel atmosphere model and lapse rate, could be in error 
by tens of mbars. This produces a systematic phase error be¬ 
tween calibrator and target and is likely the major cause of 
the poorer than expected astrometric accuracy observed dur¬ 
ing the LBC. It is expected that the addition, in late 2015, of 
more weather stations distributed over the array will improve 
the astrometric accuracy. 

5. SCIENCE VERIFICATION 

Science Verification (SV) is the process of fully testing ob¬ 
serving modes expected to be available for science observ- 
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ing by making end-to-end observations (e.g. execution of 
scheduling blocks, calibration, and imaging) of a small num¬ 
ber of selected astronomical objects. The aim is to demon¬ 
strate that ALMA is capable of producing data of the quality 
required for scientific analysis so that the observing mode can 
be offered for future science observations. To demonstrate 
ALMA’s high angular resolution capability, during the LBC 
we carried out SV observations of five targets chosen from a 
broad range of science areas (Table U). The aim was to pro¬ 
duce high fidelity, high resolution, images of continuum and 
spectral line emission using the LBC array. 

The SV targets were chosen primarily based on their suit¬ 
ability for demonstrating the long baseline capability e.g., 
having fine-scale angular structure, being less than two arc- 
sec in size, being observable at night-time during the cam¬ 
paign period, and, where possible, having previous observa¬ 
tions with other telescopes. The targets were: Juno, an asym¬ 
metric asteroid with a 7.2-hour rotation period; Mira, a well- 
studied AGB star that is the prototypical Mira variable; HL 
Tau, a young star with a circumstellar disk; 3C138, a strongly 
polarized extended quasar; and SDP.81, a high-z (z=3.042), 
gravitationally lensed, submm galaxy. Details of the targets 
and observations are given in Tableland the data are publicly 
available from the ALMA Science Porta FI Examples of the 
SV imaging results are given in three accompanying papers 
on targets HL Tau, Juno, and SDP.81 (ALMA Partnership et 
al. 2015a,b,c). Angular resolutions achieved were as fine as 
19 mas (Band 7; 344 GHz; ALMA Partnership et al. 2015b). 
In Appendix [A] we compare preliminary ALMA results on 
3C138 with a 43 GHz VLA image. Details of the imaging of 
the SV targets, including important lessons learned, are de¬ 
scribed in a CASA guide pag^3 Specific comments con¬ 
cerning the use of self-calibration to improve image quality 
are given in Appendix iBl 

6. CONCLUSIONS 

The 2014 ALMA Long Baseline Campaign achieved an 
increase of a factor of ~6 in maximum baseline length 
(~ 15 km) compared to previous test observations and a factor 
of ~10 increase compared to previous ALMA science obser¬ 
vations (a factor of ~ 100 smaller beam area). Further testing 
will be carried out in future to extend the maximum baseline 
to >15.0 km and to higher frequencies. 

Some specific results drawn from the campaign are as fol¬ 
lows. 

• Phase referencing observations should only be made 
when the short-term phase rms is < 30°, unless the 
target source is relatively compact and strong enough 
for self-calibration. This applies to all ALMA obser¬ 
vations, regardless of maximum baseline length or fre¬ 


quency. 

• Under clear skies, the WVR correction typically im¬ 
proves the phase noise by a factor of ~2. The remain¬ 
ing phase fluctuations are thought to be mostly due to 
dry atmosphere variations. 

• The prediction of short-term phase variability cannot be 
made reliably using ground-based measurements. Short 
observations of a strong source are the most reliable 
methods to determine phase conditions, as described in 
the Go/noGo procedure. 

• Systematic phase differences between calibrator and 
science target are dominated mainly by the lack of an 
accurate dry atmospheric delay model. Additional pres¬ 
sure sensors distributed across the ALMA array will in 
future improve the models. 

• The phase referencing cycle time recommended for 
long baseline observations is 60 to 90 sec between cal¬ 
ibrator observations. Shorter times do not improve sig¬ 
nificantly the image quality unless a calibrator < 1.5° 
from the target is sufficiently strong that it can be de¬ 
tected with a 6-sec integration. 

• The survey of weak calibrators will continue in order 
to increase the number of sources in the catalog and 
increase sky coverage. Alternative calibrator observing 
strategy may be needed in future in order to find the 
faintest acceptable calibrators. 

• The integration time on source may in many cases be 
driven by the time needed to obtain sufficient u-v cov¬ 
erage, rather than that needed to reach a specified rms. 
In future, detailed simulations may be needed to inves¬ 
tigate this. 

• More sophisticated methods of self-calibration may 
be needed for extended sources where the SNR on 
the longer baselines drops below that needed for self¬ 
calibration using one reference antenna. 

As a result of the extensive program of testing during the 
LBC, Science Verification at long baselines was highly suc¬ 
cessful, resulting in angular resolutions as fine as 19 mas. Ini¬ 
tial science results on the SV data are presented in ALMA 
Partnership et al. (2015a,b,c). The LBC has allowed long 
baseline (up to ~15 km) antenna configurations to be made 
available for science observations. This fulfils a major goal of 
ALMA to accurately image sources at mm and submm wave¬ 
lengths with resolutions of tens of milliarcseconds, and, to¬ 
gether with ALMA’s high sensitivity, opens up new parameter 
space for submm astronomy. 


APPENDIX 

ALMA OBSERVATIONS OF 3C138 

The source 3C138=J0521 + 1638 is a compact steep spectrum quasar with m v = 18.84 and a redshift of 0.759 (ICotton et all 
119971 ). Its angular size is about 0.4" and consists of a radio core, with a strong jet/lobe to the east and a weaker counter-lobe to 
the west. The integrated source linear polarization is 10% and its total flux density is relatively stable. 

The source 3C138 was chosen as an SV target because its angular size and small-scale structure are ideal for imaging with the 
ALMA long baselines, it is a highly polarized target, and the ALMA resolutions at Band 3 and Band 6 with a 5 to 15 km baseline 
array are comparable to that of the VLA 35 km baseline array at 43 GHz. Thus, a detailed comparison of the images made with 
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different arrays can be made. For the other LBC SV targets, the ALMA resolution and sensitivity far exceed those of other arrays 
so any detailed comparison cannot be made. Hence, the discussion here will concentrate on ALMA-VLA comparison, rather 
than any astrophysical interpretations. The analysis of the complete set of 3C138 ALMA observations (with full polarization) is 
in progress. Here, we present preliminary results FI 

The ALMA observation parameters for Bands 3 and 6 are listed in Table 2. The VLA observations at 43 GHz were made 
on February 16, 2014 in the A-configuration, and the integration time on 3C138 was 45 min. The VLA observations used 
J0530+1331 as the phase calibrator, while the ALMA observations used J0510+1800, both of which are within 4° of 3C138. The 
flux density scale for ALMA was based on the derived flux density of 1.20 Jy and 0.97 Jy (10% uncertainty) for J0510+1800 at 
97 and 241 GHz, respectively. For the VLA, the source 3C48 was used for the flux density scale. The phase referencing cycle 
time was 95 sec for ALMA and 90 sec for the VLA. 

The standard phase referencing calibration, editing, imaging, and self-calibration for the ALMA and VLA data was carried 
out using the o/?/0 software package (lCottonll2008l) . Since the structure of 3C138 is dominated by a small component, the self¬ 
calibration process was straight-forward. In order to compare the images at the three frequencies at the same resolution (91x51 
mas in P.A. -13°), each data set was weighted to include approximately the same range of spacings for each image, and then 
convolved with the above Gaussian beam size. 

The preliminary ALMA 97 & 241 GHz and VLA 43 GHz images are shown in Fig. [6] The bright, compact radio core and 
strong eastern jet and lobe respectively have spectral indices of -0.70 ± 0.03 and -0.75 ± 0.05. The western counter-jet, which 
is severely Doppler attenuated, is weak and has a spectral index of-0.95 ±0.13; its peak is just below the 3 rms intensity level 
at 241 GHz. The lowest contour level for all three images is 0.5% of their peak intensity (3 times rms), so that the peak to 
rms ratios for these images are about 500:1. The main conclusion is that the differences between the ALMA and VLA images 
are at the level of a few percent of their peak levels. The two arrays have major differences, such as their antenna, electronics, 
and correlator designs; the atmospheric conditions; and ALMA linear polarized feeds versus the VLA circular polarized feeds. 
Hence, the agreement of the images to a few percent strongly suggests that both arrays can image the radio emission from the 
sky at tens of milliarcsecond resolution with this accuracy or better. 

The ALMA Band 6 image using the high resolution data at natural weight is shown in Fig. [7] The resolution is 37x23 mas in 
P.A. -11° which is considerably higher than that used for the three frequency comparison. At this higher resolution, the western 
jet has broken into six knots and an inner jet emanating east from the core can be separated. The jet/lobe system has a slight 
curvature which is also seen on VLBA images of this source (ICotton et alJ]20Q3 ). The faint western counter-jet has a peak flux 
density of 0.25 mJy, just below the lowest contour level at 0.3% of the peak. 

SELF-CALIBRATION 

Some of the SV targets were sufficiently strong that self-calibration could be used to improve the image quality over that 
obtained with phase referencing alone. The luno images (ALMA Partnership et al. 2015a) were significantly improved with 
self-calibration and obtained a peak/image rms of typically 120 for each of the nine images, providing an increase over the 
phase-referenced only images of a factor of two to six. 

For the HL Tau continuum images (ALMA Partnership et al. 2015b), self-calibration was more challenging, because while 
the overall integrated flux is large, the source morphology is complex. Indeed, much of the disk emission is resolved out by the 
longest baselines, especially at Band 7, and for the lower frequency Bands the emission is intrinsically weaker due to the lower 
dust emissivity. Thus, the S/N for self-calibration is inadequate for the longest baseline antennas if one attempts to push to short 
enough timescales (< a few minutes) to significantly improve the phases beyond that achieved from fast-switching. Due to this 
S/N limitation on the solution interval, the self-calibration only improves the HL Tau images (peak/rms) by factors of 1.5, 1.9, 
and 1.2 at Bands 3, 6, and 7 respectively. For the much weaker source SDP.81, there is inadequate S/N to self-calibrate on a short 
enough timescale to improve the images at all (while retaining the longest baseline antennas). 

Since the 3C138 emission is dominated by a nearly unresolved core and the remaining structure is relatively simple, it showed 
the most improvement. The rms noise level decreases about a factor of 10 from the phase referenced to the self-calibrated image. 
A conservative measure is the ratio of the highest side-lobe level to the peak intensity. For the 97 GHz image, the side-lobe/peak 
intensity ratio drops from 1.4% in the phase referenced image to 0.1% in the self-calibrated image. For the 241 GHz image, the 
ratio drops from 17% to 0.6%. 

One particular complication of self-calibration at long baselines is that unless the target structure is already well-studied at high 
resolution, only a rough estimate of its correlated flux density at the longer baseline may be estimated. Therefore, in many cases 
it may be difficult to predict in advance whether a given source can be self-calibrated on the longest baselines. In future, more 
sophisticated methods of self-calibration may benefit extended sources where the SNR on the longer spacings drops below that 
needed for self-calibration using one reference antenna. Furthermore, future testing on long baselines will provide further insight 
into ALMA long-baseline imaging and self-calibration. 


This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00013.SV, ADS/JAO.ALMA#2011.0.00014.SV, 
ADS/JAO.ALMA#2011.0.00015.SV and ADS/JAO.ALMA#2011.0.00016.SV. ALMA is a partnership of ESO (representing its 
member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of 
Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The 


55 The ALMA Band 3 & 6 avergage frequencies for the initial results pre¬ 
sented here are respectively 97 and 241 GHz; only the upper sideband of the 
Band 6 data was used. 

56 Note that the ALMA data could have been processed in CASA 4.2.2 or 


higher, but was done in obit for consistency with the previously reduced VLA 
data. 
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Fig. 6.— Images of 3C138 with 91x51 mas resolution and P.A. 13° (as shown by the cross-hatched ellipse), (top) ALMA image at 241 GHz with a peak flux 
of 0.235 Jy beam -1 . (middle) ALMA image at 97 GHz with a peak flux of 0.235 Jy beam -1 . (bottom) VLA image at 43 GHz with a peak flux of 0.387 Jy beam -1 . 
For all images the lowest contour is 0.5% of the peak and the contour levels are in multiplicative increments of y/2. Details of the images are given in AppendixfAl 



Fig. 7.— Highest resolution ALMA image of 3C138 at 241 GHz. The resolution is 37x23 mas at P.A.=— 11° (shown by the cross-hatched ellipse). The contour 
levels are in multiplicative increments of y/2. The peak flux is 0.095 Jy beam -1 and the lowest contour is 0.33% of the peak. 

National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

We thank all those who have contributed to making the ALMA project possible. 

Facilities: ALMA. 
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